Direct characterization of any linear photonic device 
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We introduce an efficient method for characterizing any multi-mode linear photonic network. 
Our method employs a standard laser source and intensity measurements to directly determine all 
moduli and non-trivial phases of the matrix describing the network. We experimentally demonstrate 
our method by characterizing a 6x6 fibre-optic network and independently verify the results via 
nonclassical two-photon interference. 
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Implementation of quantum technologies requires the 
ability to realise arbitrary unitary operators, enabling ap- 
plications such as efficient quantum simulation and com- 
putation. In principle, linear photonic devices can be 
used to experimentally realise any NxN unitary opera- 
tor pQ. However, a significant remaining practical chal- 
lenge is to characterize the device once it is built. A 
known solution is to perform full quantum process to- 
mography of a device using nonclassical states [2H4] or 
coherent states [5j [6] . This standard approach however 
requires the full suite of quantum tools — such as N-mode 
quantum state preparation and measurement — and is, 
despite progress on more efficient methods such as com- 
pressive sensing [7], relatively slow and impractical for 
large interferometric devices. 

A more tractable approach, starting from the assump- 
tion of linearity, is to adapt existing methods from clas- 
sical optics. As a linear photonic circuit can always be 
cast as an interferometer with (N 2 —N)/2 beam-splitters 
PQ, it can be fully characterized by embedding it in an 
external interferometer using a local oscillator [8]. The 
interferometric stability required is sufficiently challeng- 
ing in practice that this is rarely done for quantum optics 
experiments. Recently, a method was proposed that ob- 
viates the use of an external interferometer; however, it 
requires nonclassical interference [9] for characterization 
of linear optical devices [lOj [11] . 

Here we introduce a method to characterize the NxN 
matrix representing an Af-mode linear photonic network 
using readily available standard laser source and pho- 
todetectors, and without using an external local oscil- 
lator. Our method eliminates the need for nonclassical 
interference and single-photon detectors; it is simple and 
efficient, requiring only 27V— 1 configurations for directly 
measuring all nontrivial parameters of the NxN matrix. 
We demonstrate our method by characterizing an inte- 
grated device — a 6x6 fused-fibre coupler — and highlight 
its precision by comparing measured quantum interfer- 
ence patterns with those predicted using the classically- 



estimated matrix. 

A linear multi-mode photonic network can be repre- 
sented by a matrix U that relates the creation operators 
of the input and output modes 

N 

4 = £t^4 (!) 

k=i 

For an ideal lossless network the matrix U is unitary, but 
in practice, due to loss, it is a submatrix of a unitary ma- 
trix. Notice that the knowledge of U simply determines 
the action of the network on any multi-mode coherent 
states, and in principle following the method presented 
in Ref. [6], one can predict the output state of the net- 
work for any given input state. However, for a network 
with a large number of modes, this task is believed to be 
computationally hard for certain input states [T2] . 

Generally, the elements of the matrix U are complex 
numbers Ujk=rjke ejk , where 0<r*jfc<l and O<0jfc<27r. 
Noting that the phases of the basis vectors are not phys- 
ically significant, we can absorb (2N— 1) phases into the 
basis vectors [13 . Thus, any matrix U can be decom- 
posed as a product of three matrices U=D(/j,)U D(y) , 
where D(/[x)=diag (e* Ml , e 2/i2 , . . . , e l ^ N ) and both U and 
U describe the same physical process. Without loss of 
generality, we let #ij=#ji=0, for j=l, 2, . . . , A 7 ", and we 
are left with (N— l) 2 +N 2 free parameters to be deter- 
mined in the characterization. 

Even if U is unitary, it has been shown that knowledge 
of all moduli v^k — where r 2 k is the probability of detect- 
ing a photon at the output port k when only one single 
photon was sent to the input port j — does not uniquely 
determine all phases for A">3 p~3j[14]. Therefore, in or- 
der to characterize the matrix J7, in general, we require 
probe states and measurements that are sensitive to the 
phases 6^. 

One way to achieve that, as shown in [lOj [11], is 
to insert two single photons into different input modes 
and to record nonclassical interference patterns between 
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FIG. 1: Scheme for characterizing any linear photonic device, 
U . (a) Using a 50:50 beam-splitter (BS) and phase shifter 
(0) a dual-mode coherent state, \a), is prepared and sent 
through U, where ot2—e %( ^ ot\. By sequentially inputting \a.2) 
into modes 2,3,...,iV, and varying the phase over at least 2tv, 
all phases of matrix U can be directly measured using photo- 
diodes (PD). (b) Experimental realisation. The device-under- 
test is a 3x3-mode fused- fibre beam-splitter (FBS), which 
constitutes a 6x6 photonic network when polarization is in- 
cluded. Orthogonal polarisation modes are resolved using fi- 
bre polarization beam-splitters (FPBS) at its outputs. Inter- 
ferometric probe states between pair- wise input combinations 
are prepared with two polarization beam displacers (BD), and 
half- wave plates (HWP). Our bulk optics method to direct co- 
herent states here is used for convenience only. Finally, the 
outputs are monitored with fast photo-diodes connected to 
an oscilloscope (OSC) while the phase (j) is scanned with a 
motorised linear micro-translation stage. 



different combinations of output modes. Alternatively, 
the nonclassical interference can be simulated using two- 
mode coherent states with randomized relative phases, 
but at the cost of an increased level of noise and reduced 
classical interference visibility [TSJIIS]. Our method takes 
a more direct route, only requiring a standard laser 
source, split at a beam-splitter, with a varying relative 
phase between the resulting dual-mode coherent state; 
see the experimental setup in Fig. [IJa). An interesting 
feature of this method is that it enables us to measure 
all the nontrivial phases of the matrix U directly without 
solving any trigonometric equations. It thus significantly 
simplifies the task of characterization. 

Our method works as follows (details of the algorithm 
can be found in the Appendix): 



1. Send coherent state of light into the N input modes 
individually, and measure the output r 2 - k of all 
fc=l, 2, . . . , N outputs modes for every input j. 
From this, obtain r^. 
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2. Prepare a dual-mode coherent state into pairwise 
input combinations between input 1 and the input 
j=2, 3, . . . , A/", with a variable phase (j). Measure 



FIG. 2: Experimental characterization of a linear photonic 
device, (a) The moduli r 2 jk of the experimentally measured 
Uexp- The x and y axes correspond to the input and outputs 
modes, j and k respectively, (b) Representative experimental 
data for obtaining Ojk- Injecting coherent state of light into 
input modes 1 and 3 of the 6-mode device. The amplitudes 
(voltage at output photo-diodes) of the six output modes (1-6 
from top-to-bottom) oscillate as the phase (j) is swept in time. 
The red and blue lines are measured data and theoretical 
fits to Acos(0 — Ojk) respectively, (c) All phases Ojk of the 
measured matrix [/ exp . The entire characterization method 
was performed 10 times to obtain experimental uncertainties, 
however in all plots the experimental uncertainties are not 
visible on the scales shown. 
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FIG. 3: Independent verification of the measured matrix U exp . (a) Experimental schematic. A pair of photons is generated via 
spontaneous parametric downconversion (SPDC) in a nonlinear beta-barium-borate (BBO) crystal pumped with a mode-locked 
pulsed laser at 410 nm with type-I phase matching to produce polarization degenerate photons at 820 nm. After being spectrally 
filtered (FWHM 2 nm) the individual downconverted photons are steered into the optical modes of the linear photonic network 
by a series of beam displacers (BD) and half- wave plates (HWP). Upon passing through the photonic circuit the photons are 
detected using avalanche photo-diodes (APD) who's coincident signals are monitored using a commercially available counting 
logic, (b) Measured nonclassical visibilities vs. predicted visibilities are shown in the left and right panels are for photons input 
into modes {1,2}, and {1,6} respectively. The red bars show the directly measured nonclassical visibilities, the blue bars show 
the predictions from the measured matrix [/ exp , and errors are given at the top of each data point. The numbers on the x-axis 
show the corresponding output modes for which that nonclassical interference is measured. 



the output intensities of all modes k=l, 2, . . . , N as 
a function of <fi. This yields all non-trivial phases 
Ojk- 

We tested our method by characterizing a fused fibre- 
optic beam-splitter, the experimental setup is shown in 
Fig. [ljb). The network itself is composed of one 3x3 
non-polarizing fused fibre-optic beam-splitter (FBS) with 
three 2x2 orthogonally-polarizing beam-splitters at each 
of its output modes. By mapping onto orthogonal polar- 
izations at the input of the initial FBS the whole network 
is described by a 6x6 matrix Z7 exp . The input modes 
are labelled {1, 2, . . . , 6}={|#>i , \V) 1 , \H) 2 , . . . , \V) 3 }, 
where \H) 1 is the horizontally polarized mode for spa- 
tial mode T of the FBS. 

In the setup we used a series of polarization beam 
displacers and half-wave plates to prepare input probe 
states, allowing for phase-stable interferometric measure- 
ments and polarization control for the input of the FBS. 
The phase (j) was controlled by a motorised linear micro- 
translation stage at input mode '1' to introduce an opti- 
cal path difference of 0.1 mm, at a speed of 0.05 mm/s, 
between two inputs. Scanning over this short time win- 
dow limits the effect of thermal drift on the classical 
interferometer, therefore removing the need for active 
locking. The outputs were coupled to fast photodiodes 
and monitored simultaneously on an oscilloscope while 
the phase <ft was scanned. All characterization measure- 
ments were performed with a 100 jaW laser diode filtered 
to have a centre wavelength of 820 nm and a full-width 
half-maximum bandwidth of 2 nm. 

We first measured the 36 output intensities for the six 
individual inputs shown graphically in Fig.[2|a). We then 
recorded interference fringes for the pair- wise input com- 
binations discussed above, and fitted sinusoidal curves to 



the resulting photocurrents to obtain experimental val- 
ues for Ojk (see Appendix for details), see Fig|2|b) and 
(c). From these measurements, we reconstructed the 6x6 
matrix U exp . 

We verify our experimentally obtained matrix [7 exp by 
measuring two-photon interference inside the linear op- 
tical network [9j. As this is a fourth-order interference 
effect, it provides a suitable independent verification for 
the validity of U exp which we obtained with our second- 
order interference method. We created and sent two 
single-photons into the 6x6 network, and measured the 
coincidences at all fifteen (6 choose 2) pairwise combi- 
nations of output modes using the experimental setup 
shown in Fig. [3^a). The nonclassical visibility is calcu- 
lated as (C ou t-C in )/C ou t where C in and C out are the 
coincidence count rates of photon pairs at the outputs of 
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FIG. 4: The matrix V exp Vj xp . The diagonal hatched squares 
are equal to unity by construction (see Appendix) and the off- 
diagonal elements are coloured according to their value given 
by the color bar on the right. 
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the circuit inside (maximum temporal overlap between 
photons) and outside (no overlap) the interference dip, 
respectively. The results for two different input configu- 
rations are shown in Fig.[3jb). The obtained interference 
patterns are in excellent agreement with those predicted 
by Z7 exp . 

Having verified our method, we will now discuss how 
U ex p — which we do not expect to be a unitary matrix be- 
cause of inevitable loss — can be embedded into a larger 
matrix which will more closely approximate a unitary 
matrix. Assuming the optical loss is equal for different 
paths connecting specific inputs to outputs, we can model 
this loss as virtual beam-splitters on the input modes. 
With this we obtain the 12x12 matrix V exp following 
the algorithm described in the Appendix. If our linear 
photonic device satisfied these assumptions, we would 
expect V exp to be unitary. A simple measure of this uni- 
tarity is to check how close V exp V^ p , given in Fig. [4| is 
to identity. In our case KxpKxp nas non-zero but very 
small off-diagonal elements. This is indicative of both 
inevitable experimental uncertainty in the measured ma- 
trix elements and the effect of path- dependent loss in the 
linear optical network. 

The question remains as to what unitary best describes 
the linear optical network. It has been shown that the 
closest unitary matrix to V^ xp can be found by using the 
polar decomposition t/=(V r e xpKxp) _ ^^exp [17 . The re- 
sulting unitary matrix U does not noticeably alter the 
predicted two-photon interference visibilities of Fig. [3jb) , 
and therefore describes our device with a good approxi- 
mation. 

As photonic quantum technologies mature beyond 
small-scale proof-of-principle demonstrations [T8H2Q] . 
there is increasing requirement for improved methods 
for full tomography, process validation and verification. 
Areas of direct applicability include the experimental 
characterization of large linear waveguide arrays for pho- 
tonic quantum walks [2TJ [22] , especially those built in 3D 
where current approaches using top-down imaging are 
not possible [23]. One of the most exciting applications 
for larger systems pertains to an intermediate model of 
quantum computation — namely, that of BosonSam- 
pling [12]. Our method provides an efficient means for 
characterizing large linear optical networks to obtain 
the scattering probabilities of multi-photon processes, 
which is a crucial component of a photonic BosonSam- 
pling experiment. For these larger characterizations a 
many-mode optical switch board built into inherently 
stable integrated wave-guide circuits [24] could be used 
to prepare and steer the required dual-mode coherent 
states, as opposed to a larger bulk-optics setup. A 
simpler way to achieve the required small displacements 
is to tilt optical components in the probe beam path [16] . 
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Appendix 



Here we show how to construct the matrix V exp (la- 
belled V here for neatness) directly from the experimen- 
tal measurements. We model a lossy NxN linear pho- 
tonic device by adding virtual beam-splitters with trans- 
missivities rjj (0 < rjj < 1) to each input of an NxN 
unitary matrix with elements Ujk=Tjk^ %ejk ^ and the re- 
sulting 2Nx2N unitary matrix is given by 



V 



fjU 7] 



(A.2) 



where 



r/ = diag(77i,r/ 2 ,...,w) 5 



V = diag Jl - rjf, a/1 - r/|, . . . , Jl - rj 2 N 



Using Eq. ([I]), the output amplitudes of a general N- 
mode input coherent state \ql\, ol2-> • • • ? &n) are given by 



N 



(3k = ^2v jk aj. 



(A.3) 



We can thus send a coherent state with known intensity 
I to input port j, where all other input modes are in 
vacuum state, and measure the intensity from all output 
ports simultaneously 

Ik=HVjr jk )\ fc = l,2,..., N. (A.4) 
We also have 



1 N 
(%) 2 = yE^ 



(A.5) 



k=l 



Thus using equations (A.4) and (A.5) for j=l, 2, . . . , N 



all the moduli can be obtained. 

The next step is to measure all nontrivial phases using 
two- mode coherent states. We send a coherent state \a) 
to a 50:50 beam-splitter and use a phase shifter to control 
the relative phase between the output states \a±) and 
|a2) with the same intensity /; see Fig. [TJa). We send 
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|ai) to input mode T and e 1 ^ \a2) to input port j. The 
intensities of the output coherent states are given by 



I \Vik+V jk e 



icj) I 



(A.6) 



As all elements in the first row and the first column are 
real, the above equation becomes, for k=f=\ 



h = l 
= I 



m r ik +Vj r jke 



{r}ir lk f+ (Vjrjkf+ZViVjrikrjk cos( 



9jk) , 
(A.7) 



and for k=l 

h = i 



(m r ii f + (Vjrji) 2 + 2mVj r ii r ji cos(4>) 



(A.8) 

When I\ attains its maximum value we have 0=0, and 
(j)=7r for its minimum value. This serves as our reference 
mode, and without loss of generality we always choose 
I\ at its maximum for our reference of <p=0. Having 
this capability, we further sweep (j) until I k attains its 
maximum value (see Fig. [2jb)) and using Eq. (A.7) the 
unknown phases can be found as 



2tt- 



(A.9) 



Therefore, by repeating this procedure for \a.2) input into 
mode j=2,3, ...,N all the nontrivial phases 6j k of the 
matrix V can be directly measured. 
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